Likelihood Ratio Introduction 


A new feature added to the Generic Monte Carlo Engine with Volatility Skew is the Likelihood 
Ratio method for calculating risk sensitivities. The Monte Carlo engine is used for pricing and risk 


management of various products 


The Local Spline Dupire volatility calibration serves its intended purpose to produce a local 
volatility surface which can be used to calculate the price and risk sensitivities of various products, 
based on Monte Carlo simulations or the PDE approach. The Likelihood Ratio method serves its 


intended purpose to calculate risk sensitivities in the Monte Carlo approach. 


The key model assumptions are: 
a. The asset values can be accurately expressed using a volatility skew model. 
b. The interest rates are deterministic. 
c. The volatilities of the exchange rates are flat (constant). 
d. The pairwise correlations between assets are constant. 
e. The pairwise correlations between assets and exchange rates are constant. 


f. The future dividend dates, amounts and yields are known and fixed. Dividend yields are 


treated to preserve the dividend cash-flows. 


The main risk uncertainties are the implied volatility surfaces of the underlying assets. Other risk 


uncertainties are correlations among assets and between assets and exchange rates. 


The Local Spline Dupire volatility calibration can be used to price various products based on either 
Monte Carlo simulations or the PDE approach. The current implementation of the Likelihood Ratio 
method can be used when the parameter IsDistribution is set to True and the payoff function 


depends on the asset prices only at maturity. 


The accuracy of the price and risk sensitivities produced by the Monte Carlo engine depends on the 
number of simulations. The convergence levels for each specific product will be estimated 


separately. 


The generic Monte Carlo engine with volatility skew generates random paths for one or more 
underlying assets. These paths are used to calculate the price and risk sensitivities of various 
contingent claims based on these underliers. The model incorporates the volatility skew by 
assuming that the volatility is given either by a function of time to maturity and spot (local 


volatility) or by an implied distribution at maturity, depending on the simulation method. 


The generic Monte Carlo simulation produces a set of random paths which are assumed to be 


representative of the probability space defined by the full range of possible asset prices (see 


https://finpricing.com/FinPricing-ProductBrochure.pdf). 


The Local Spline (LS) calibration assumes that, in the log-moneyness space, the local volatility 
dependency on time is piecewise-constant, while the dependency on strike is given by a Bezier 
curve. Thus, local volatility is defined as a set of curves, one for each maturity in the implied 
volatility table. The curves are defined by uniformly spaced control points between the lowest and 


highest available strikes. 


Users can specify the number of control points and a penalty term, P = 0, which restricts the arc 
length of the curve in order to reduce oscillations and enforce smoothness. Each curve is the 
solution of a least-squares optimization problem that minimizes the differences between model call 


prices and market call prices. The model call prices are determined by solving a numerical PDE for 


either the Fokker-Planck equation for the probability density function or the Dupire equation for 


call prices. 


Under the sticky strike rule, when the current spot differs from the volatility spot, the surface is 
recalibrated. This introduces numerical instability since the minimizing routine may reach a 


different solution. 


The new method, named LS Dupire, avoids the recalibration step by deriving the local volatilities 
from the PDE call grid, using the Dupire formula. The Dupire equation in the log-moneyness space 


is ([5]): 
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where 
o= o(k, t) the local volatility 
k= el | the strike in the log-moneyness space 
Ot 
K the strike 
Fo, the forward price from time 0 to time ¢ 
c=C (k ot ) = UK, ‘ the normalized call price 
Fo, 
C (K i t) the undiscounted Black-Scholes call price 


For the new method, the volatility is calibrated to the volatility spot using the same objective 
function as the LS calibration. The new method can fit the volatilities using either Bezier curves or 
monotone cubic splines (the latter being preferred). The control points for the monotone cubic 


splines are the available market quotes (which are not uniformly spaced). If the current spot differs 


from the volatility spot, implied volatilities are derived from the call prices on the PDE mesh from 
the initial calibration. The local volatilities for the current spot are obtained from the implied 


volatilities using the Dupire equation: 
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where 

= o(K : t) the local volatility 
v=v(K,t) the implied volatility 
K the strike 
r(t) the interest rate 
q(t) the dividend yield 


. ae term coming from the Black Scholes formula 
S the current spot price 
The new method is stable when the market changes because the spot affects only the d term. 


Moreover, the resulting local volatility surface is smooth because the implied volatilities are 


arbitrage-free. 


The Monte Carlo engine calculates risk sensitivities using the finite difference (FD) approach. This 


method is simple to understand and implement but it has a few well-known drawbacks: 


e It requires a large number of simulations to achieve reasonable convergence. 


e Each risk sensitivity requires at least one additional valuation. 
Thus the FD approach is slow and computationally demanding. 


The Likelihood Ratio method (LRM) assumes that the price depends on the selected risk factor 
only through the probability distribution. More precisely, consider the value of an option on a 


basket of n assets: 
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where fi is the payoff function which depends only on the spot prices at maturity, and f is the 
terminal joint distribution of the assets S,,S,,...,5,. Let 0 be a variable that f depends on but 


h does not. Then, the price sensitivity with respect to @ is: 
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Thus, the risk sensitivities are estimated along the original simulation paths using a payoff function 


Clog f 


adjusted by the likelihood ratio . This makes LRM more efficient than FD and the resulting 


sensitivities more stable to the number of simulation paths used. The Appendix provides further 
details regarding the LRM for the Gaussian copula as well as comparison tests between LRM and 


FD. 


It should be mentioned that the current implementation of LRM can be used only when the 
IsDistribution parameter is set to True and the payoff function depends on the spot prices only at 


maturity. 


The Local Spline Dupire volatility calibration can be used to price various products based on either 
Monte Carlo simulations or the PDE approach. The current implementation of the Likelihood Ratio 
method can be used when the parameter IsDistribution is set to True and the payoff function 


depends on the asset prices only at maturity. 


The graphs below show the changes in PDE price and risk sensitivities with respect to the spot 
value. The changes in price and discount rho are minimal, while changes in gamma and vega are 


substantial: 
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PDE Delta 
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PDE Gamma 
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PDE Vega 
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Recall that under LRM the risk sensitivities are obtained by adjusting the payoff function: 


Similarly, second order derivatives are given by: 


av a? 
= A\S,,S,,...55 ] S305 5en5 0. S,,S85,...59, JdS 
06,08, J ( —— Le 08 f( eee dye iors ) 3 


00, 


é é 
+ | HS Sre08,|{ ts AS S308) { SBI S205) PES 
1 


where f, (s ) are the marginal densities, F, (s ) the marginal cumulative distributions and 


c(F.(S, ; F(S, ), eae S )) the copula term. The Gaussian copula is: 


n 


where & is the correlation matrix and u, = N a (S ; ). 


In this case, the likelihood ratio becomes: 
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When the parameter IsDiagonalGamma is False, gamma is calculated as shown below. Let 


f(S,,a)=f(S, TOD asa deO feasdGO +ab,), where b,, j #i represents the change in S, with 


respect toa move in S;,. Also let: 


far = F(S,(1+e)A) 
fi = F(S; 1) 
fin =f (S,-1) 
faa = f(S;(1-€)-1). 
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The likelihood ratio for the gamma we want is log f + [2 oe ie {2 log f ) which is 
i a 


estimated as follows: 


bre log( fi, )—log(f;, )—log(Fj, )+ log (fi) 
aS,da («s,) 


3 alos (r;,)-log(fj, )+log(Fj, )-log(F,,)) 
as, &S, 


a ; (log(s,, )+log(F;, )-log(F;, )-log( i, )) 
0a &S, 


